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This  report  is  a  study  of  the  relative  accuracy  of  six  alternative  methods 
of  calculating  a  solution  involving  shocks  and  a  rarefaction  wave  for  the 
equation 

ut  +  (•j  °2  +  *0  =  0  • 

The  equation  is  equivalent  to  that  proposed  by  Burgers  as  a  simplified  model 
for  shock  problems  in  fluid  dynamics.  The  results  are  believed  to  be 
suggestive  as  to  calculations  in  such  problems. 
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*  A  COMPARISON  OF  NUMERICAL  SCHEMES  TO  CALCULATE  THE  SOLUTIONS 

OF  A  NON-LINEAR  PARTIAL  DIFFERENTIAL  EQUATION  WITH  SHOCKS 

INTRODUCTION 

A  simple  problem  with  shocks  and  rarefactions  for  the  equation 

(1)  ut  +  (iu2  +  u)x  =  0 

has  been  subjected  to  calculation  by  various  procedures  and  all  the  results 
compared  with  the  exact  solution.  The  methods  thus  tested  are 

(1)  Lax's  original  centered  difference  method  [l], 

(2)  a  left  difference  method  without  viscosity  [2], 

(3)  the  viscosity  method  of  Lax  and  Wendroff  [3] 

(3.0)  without  added  artificial  viscosity, 

(3.1)  with  added  artificial  viscosity, 

(4)  a  modification  of  the  Lax-Wendroff  method  for  which  convergence 
has  been  proved  [4] 

(4.0)  without  added  artificial  viscosity, 

(4.1)  with  added  artificial  viscosity. 

The  programming  of  these  schemes  for  an  I  EM  704  machine  was  performed  mainly 
by  W.  Parr,  with  the  assistance  of  Mrs.  S.  Madigosky,  for  whose  patience  and 
intelligent  care  I  wish  to  express  my  gratitude  and  thanks. 

Only  calculations  of  roughly  the  same  cost  are  compared.  A  multiplica¬ 
tion  or  division  costing  about  the  same  as  ten  additions,  the  measure  of 
relative  cost  has  been  taken  as 

A  C  -  +  a 

10,000  hk  ’ 
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where  h  and  k  ara  tha  horisontal  and  tha  vartieal  distances,  respectively, 
between  consecutive  grid  points  (thus,  l/hk  measures  the  total  number  of 
calculations  performed  in  unit  area  of  the  x,t-plane),  p  denotes  the  number 
of  multiplications  and  divisions  performed  at  &  typical  grid  point,  and  a 
the  number  of  other  operations  at  that  point. 

The  problem  treated  is  that  of  finding  a  generalised  solution  ^  of  (l) 
with  initial  data  prescribed  as 


u(x,0)  = 


0  for  x  <  0 
2  for  0  <  x  <  0.3 
1  for  0.3  <  x  <  0.9 
0  for  x  >  .9 


This  problem  is  easily  solved  explicitly  with  a  rarefaction  wave  centered 
at  (0,0)  and  shocks  issuing  from  (.3,0)  and  (.9,0),  each  shock  being  a 
line  of  slope 


u.  +  u 


=  1  + 


where  u+  and  u_  are  the  limiting  values  of  the  solution  at  the  right  and 
at  the  left  of  the  shock,  respectively.  In  the  interval  0  <  t  <  0.6,  in 
particular,  the  solution  of  the  problem  is  given  by 

/  0  for  x  <  t 


u(x,t)  = 


^  -1  for  t  <  x  <  3t 

2  for  3t  <  x  <  .3  +  |  t 

1  for  .3  +  |t  <  x  <  .9  +  |t 

0  for  x  >  .9  +  |t 


(See  Figure  1.)  This  solution  is  known  to  be  unique 
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Note  that,  in  this  problem,  0  <  u  <  M,  where 
M  =  2  . 

The  same  bounds,  by  a  priori  reasoning,  applies  to  the  calculated  values  of 
u  in  the  first,  second,  and  fourth  schemes,  and  for  present  purposes  will  be 
assumed  in  the  third.  Note  also  that  all  schemes  considered,  ex&ept  the 
third,  have  been  proved  to  render  approximations  to  the  solution  which  are 
as  close  as  desired,  provided  h  and  k  are  sufficiently  small  and  the  ratio 
■0  =  k/h  appropriately  bounded. 

THE  CALCULATION  SCHEMES  TESTED 

The  notation  in  the  following  is  that  of  [4]:  v^  denotes  the  value  of 
u  ,  as  calculated  according  to  the  scheme  in  question,  at  the  grid  point 
X  =  lh,  t  =  Jk,  =  ( Vlij  -  r1})A  and  Wjj  =  h.^  =  vJ+1>J  -  ( 

F(u)  =  u2/2  +  u  ;  =  F(v^j ) ,  f'aj  =  F'tk^j)  .  Also,  as  above,  $=  k/h  . 

1.  Lax's  original  centered  difference  method.  Lax's  original  scheme 
was  based  on  the  difference  equations 

Vi,J+l  =  2(vi+l,j  +  Ti-l,j*  ”  2(Fi+l,j  “  Pi-l,J^  * 

According  to  N.  D.  Vvedenskaya  [6]  (Theorem  1),  the  scheme  converges  when 
•0  max|F'(v)  |  <  1  ,  or,  in  our  case,  when  -0  <  l/(l  +  M)  =  l/3  . 


2.  A  left  difference  method  without  viscosity.  In  this  scheme. 


ij+l 


=  v 


ij 


-  «(F, 


ij 


"  Pi-l,J} 


i 


m 
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The  asymmetry  of  the  scheme  is  permitted  because  of  the  non-negativity  of  F' . 
According  to  a  remark  at  the  end  of  the  appendix  below,  convergence  occurs 
when  -0  <  1/4  . 


3.  The  viscosity  method  of  Lax  and  Wendroff.  In  this  method, 

2 


v  =  v  -  -(F  -F  )  +  -f  (F1 

i,j+l  Vij  2  l+l,  j  *i-l,j;  4U  ij 


Fi+l,j)(vi+l,j  "  vij)  " 


wnere 


Qij  =  “  F. 


ij 


(v 


i+l,j'  v*i+l,j  ’ij 


-  O  »  B 


I 

being 


a  non-negative  constant. here  taken  as  zero  in  method  3-0  and  as  1/4  in 
method  3.1.  Lax  and  Wendroff  [3,  p.  227]  expect  the  scheme  should  be  stable, 
in  the  case  B  =  l/4  ,  when  -0  max|F|j|  <  .78  ,  thus  when  -0  <  ,78/(l+M)  = 
.26.  Stability  in  the  case  B  =  0  by  their  formula  demands  merely  -0  <  1/3  . 
(The  B  in  this  paper  is  one  fourth  the  quantity  Lax  and  Wendroff  call  B.) 


4.  A  modified  Lax-Wendroff  scheme.  The  calculation  scheme  here  con¬ 
sidered  is 


'i,j+l~vij  ^Fij“Fi-l,j)+  2  ^Fij2(vi+l,j"Vij)"Fi-l,j(Vij~Vl,j)^+e(Qi-l,j  Qij): 


the  constant  B  entering  Q,^  being  taken  as  zero  in  Method  (4.0)  and  as 
1/4  in  Method  (4.1).  For  B  >  0,  it  is  proved  in  [4]  that  this  method 
converges  if  -0  is  sufficiently  small,  but  larger  possible  choices  for  ■O 
than  emerge  from  that  discussion  are  determined  in  the  appendix  below;  these 
are 
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«  =  .16  for  B  =  0  (Method  (4.0)) 

=  .156  for  B  =  1/4  (Method  (4.1)). 

Two  cost  levels  (C  =  8  and  C  =  32)  were  arbitrarily  fixed  for  all 
the  above  schemes  and  values  of  4  selected  near  their  permitted  maxima. 
Then  mesh  widths  h  and  k  were  obtained  to  accord  in  each  case,  with  the 
predetermined  C  and  -6  .  The  values  of  the  parameters  thus  entering  are 
displayed  in  Table  1. 
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CONCLUSIONS 

The  results  of  the  calculations  by  each  aethod  except  4.1  were  tabulated 
after  some  thinning  for  the  rectangle 

0  <  x  <  1.8 
Rs  0  <  t  <  .6  , 

the  case  C  =  32  being  presented  here  in  Tables  4-8.  The  Bean  absolute 
error  and  root  mean  square  error  over  all  mesh  points  in  R  and,  separately, 
in  its  top  boundary  line  L  are  given  for  all  calculations  in  Table  2.  These 
results  suggest  the  following  observations,  at  least  with  respect  to  the 
example  studied: 

1.  The  use  of  artificial  quadratic  viscosity  (B^O)  in  Method  3.1  did 
not  clearly  reduce  the  average  error  and  in  Method  4.1  greatly  increased  it. 
Other  values  of  B  and  lower  mesh  width  ratios  might,  however,  have  led  to 
better  results. 

2.  Lax— Wendroff  viscosity  worked  well  in  the  method  of  centered  differ¬ 
ences  (3.0)  for  which  it  had  been  proposed,  but  disappointingly  in  the  leftr- 
difference  scheme  (4.0).  Perhaps  the  latter  aethod  would  have  performed 
better,  however,  at  a  higher  cost  level  at  which  the  x-axis  would  be  more 
finely  subdivided. 

3.  Quadrupling  the  cost  less  than  halved  the  average  error. 

4.  Methods  2,  3-0,  and  3.1  worked  best:  at  identical  costs,  closely 
comparable  mean  errors  resulted  from  the  simple  scheme  and  refined  mesh  of 
Method  2  and  the  refined  schemes  and  coarser  mesh  of  Methods  3.0  aad  3.1. 

Additional  conclusions  arise  from  comparing  the  three  most  successful 
methods  on  a  region  of  relatively  great  transitions. 
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C  =  8.00 


AVERAGE  ERRORS 


C  =  32.00 


*R 

*L 

SR 

SL 

.20 

.18 

.31 

.28 

.14 

.12 

.23 

.20 

.16 

.11 

.29 

.22 

.18 

.11 

.29 

.18 

.24 

.21 

.34 

.32 

*r 

■l 

SL 

.12 

.11 

.20 

.18 

.08 

.08 

.15 

.16 

.08 

.07 

.18 

.23 

.09 

.08 

.18 

.19 

.16 

.15 

.25 

.25 

.34 

.34 

.55 

.63 

:  x  <  i 

.8,  0 

<  t 

<  .6. 

0  <  x 

<  1.8, 

t  = 

:  [.6]k,  where 

.  multiple  of 

k 

that  does  not 

exceed  a  . 

Ep  =  average  absolute  deviation  of  calculated  fro*  exact  values 
of  the  solution  at  the  lattice  points  in  R  . 

=  average  absolute  deviation  of  calculated  fro*  exact  values 
of  solution  at  the  lattice  points  on  L  . 

3p  =  root  mean  square  deviation  at  the  lattice  points  in  R  . 

SL  =  root  mean  square  deviation  at  the  lattice  points  on  L  . 

C  =  index  of  total  cost  . 


Better  results  were  obtained  by  Method  4.1  with  B  =  1  and  a  value 

of  9  much  <  -0  ;  they  still  compared  poorly  with  the  others, 

wr 


Table  2 
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For  C  =  32  ,  in'Tables  9,  10,  and  11  we  present  with  no  thinning  in  the 
x-direction  such  a  comparison  over  the  trapezoid 

-.05  +  3t  <  x  <  .95  +  3t/2 
Ts 

.3  <  t  <  .6  , 

this  area  covering  much  of  the  vicinity  of  the  two  shocks  and  a  small  part 
of  the  rarefaction  wave.  The  mean  absolute  error  and  the  root  mean  square 
error  over  all  the  mesh  points  of  T  appear  in  Table  3*  From  these  compila¬ 
tions  it  seems  reasonable  to  add  the  following  conclusions: 

5.  By  all  three  methods  (2,3*0  and  3*1),  calculated  values  tend  to  be 
too  low  to  the  left  and  too  high  to  the  right  of  a  shock.  These  methods  thus 
dull  the  apparent  sharpness  of  a  shock,  but  Method  2  perhaps  to  the  greatest 
extent. 

6.  Method  3*0  seems  most  prone  to  occasional  wild  errors  (occurring  in 
this  example  to  the  right  of  the  first  shock) ;  these  errors,  however,  are 
rapidly  corrected. 

7.  The  average  absolute  error  over  T  is  about  the  same  for  the  three 
methods;  the  mean  square  error,  however,  is  appreciably  less  under  Methods  2 
and  3*1  than  under  Method  3.0. 

Further  study  of  our  six  methods  might  consider  the  effects,  at  constant 
cost,  of  varying  the  mesh  width  ratios  below  their  maximum  permissible  values 
and  also  of  varying  the  coefficient  of  artificial  viscosity  B  . 
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AVERAGE  CALCULATION  ERROR 
IN 

REGION  OF  RAPID  TRANSITIONS 


Method 

C  = 

32 

et 

ST 

2 

.25 

.29 

3.0 

.23 

.36 

3.1 

.23 

J 

.29 

T  denotes  the  trapeaoid  -.05  +  3t  <  x  <  .95  +  3t/2 

.3  <  t  <  .6 

Et  =  average  absolute  deviation  of  calculated  fro*  exact  values 
of  the  solution  at  the  lattice  points  in  T  . 

ST  *  root  *ean  square  deviation  at  the  lattice  points  in  T  . 

C  a  index  of  total  cost 


Table  3 
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TABLE  4  CENTERED  DIFFERENCE  METHOD  OF  LAX!  METHOD  I,  C  =  32. 
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TABLE  8  MODIFIED  LAX- WENDROFF  METHOD  WITHOUT  ADDED  VISCOSITY*.  METHOD  4.0  C*32. 
CALCULATED  VALUES 
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TABLE  9  LEFT  DIFFERENCE 
EXACT  MINUS  CA 


*  The  dashed  line  is  the  right  hand  boundary  of  the  rarefaction  wave:  the  solid  lines 
intersecting  it  are  the  two  shocksj.  The  coordinates 

x  =  .02739  (1  -  1),  t  =  .066845  (J  -  1) 


to  which  any  particular  entry  corresponds  are  ascribed  to  the  decimal  point  for  that 
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(OD  WITHOUT  VISCOSITY:  METHOD  2,  C- 32. 
TED  VALUES  IN  REGION  OF  RAPID  TRANSITION* 


TABLE  II  LAX-WENDROEF  METHOD  WITH  ADDED  VISCOSITY:  METHOD  3.1,  C  =  32,  B  =  l/4 
EXACT  MINUS  CALCULATED  VALUES  IN  REGION  OF  RAPID  TRANSITION* 
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(1-1),  t  =  .01258  U-l) 
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APPENDIX 

Our  aim  here  is  to  find  as  large  values  of  -0  as  possible  for  which 
methods  4.0  and  4*1  will  converge.  Two  procedures  will  be  compared,  the 
first  that  described  in  [4]  and  the  second  a  modification  of  it.  The  second 
procedure  turns  out  to  be  the  better  and  leads  to  the  values  actu. ly  used 
above. 

Calculation  scheme  4  is  that  of  [4]  with 

f=r,  g  -  0,  .  =  0,  ^  «  -f(i  +  *  Qu  > 

3/ 

see  pp.  5  and  21  and  equation  (4.4b)-ri-/  In  inequality  (4.7),  in  particular, 
therefore,  as  we  see  at  once, 


bij  =  +  Vi-l,j)2  +  BlWi-l,j  *  »  cij 


=  |(1  +  vij)2+  B  jW 


ij 


dij  =  eij  =  0  ' 

To  obtain  inequality  (4.8),  it  is  well  to  note  that 

Qij  =  “  BlWijJWij  =  f(Wij)  ’  Wher®  f(W)  =  "  BIWIW  » 

and  to  use  the  method  described  for  case  (4. 5c).  In  the  result,  we  evidently 
will  have 

Bu  =  !(1  +  vij)2*  *iMij  1  •  °ij  =  !(1  +  vi,j>2  +  2Biwij  1  • » *  “  hi  ■ 


Hu  = 1  +  (,ij +  -  Du  *  ■  Eij  *  Eij  ■ 0  - 


where  W ! .  and  W" .  are  numbers  between  V. .  and  W.  ,  and  between  U. .  and 
ij  ij  ij  i-l,j  ij 


W 


i+l,j,  respectively. 
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Both  ways  of  bounding  ■©  begin  with  the  argument  for  Theorem  5.1.  In 
this  argument,  <6  is  restricted  by  the  single  requirement  that  the  co¬ 
efficient  of  v^j  in  the  right  member  of  (5.2)  be  non-negative,  and,  there¬ 
fore,  that 

1  -  -0(F'*  +  b*  +  c*  =  ks«)  >  0  , 

the  star  in  each  case  indicating  a  maximum  value  of  the  quantity  concerned 

for  the  various  possible  values  of  its  indices  or  arguments.  Since 

F,#  =  1  +  M,  b*  =  c*  =  ^(1  +  M)2+  MB  ,  s*  =  0,  we  have  the  condition 

(2)  (1  +  M)2-©2  +  (1  +  M  +  2MB )£  -  1  <  0  . 

M  being  equal  to  2,  for  B  =  0  we  thus  have 

(3) c  «  <  (51//2-l)/6  =  .206+ 

and  for  B  =  l/4 

(3)i  «  <  =  -178- 

In  the  first  of  the  two  ways  of  bounding,  -0  now  is  subjected  to  an 
additional  restriction  arising  from  the  argument  for  Theorem  5.2.  According 
to  this  argument,  the  coefficients 

=  1  -  «i+v  !<1+Ti,j)2t  !(i*vi,j)2+  2BiwiJi +  i 

and 

a  =  |(1  -  20(1+M)) 

in  (5.8)  (b=0  in  this  equation)  must  be  made  to  satisfy  the  conditions 
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(p.14)  p^j  >  2/3  and  a  >  0  5  hence,  we  require 

(1  +  M)202  +  (1  +  M  +  4BM)«  -  1/3  <  0  and  -0  <  l/2(l+M)  . 

The  first  of  these  inequalities  would  imply  the  second.  Since  M  =  2,  the 
first  inequality  is  equivalent  to 


(4)a 


0 


.  (7/3)l/2  -  1  _ 

<  -y  — —  —  • 


087  for  B  =  0 


and  to 

(4)1  ^3—  =  .06  for  B  =  1/4 

It  expresses  the  first  of  the  two  alternative  bounds  we  have  considered  upon 
0  . 

The  second  method  Of  bounding  0  comes  from  a  second  method  of  proof 
of  Theorem  5.2  in  which  an  argument  of  Vvedenskaya  has  been  applied.  We 
begin  with  inequality  (5.7)  which,  in  the  present  case,  reads 


<  t1  -  e(Fij +  Bu +  +  «<Fij +  Vvi.j +  ecijVi,j 


V,  >v 


_  —I  f*  _  20(1  +  2  k  pn  2 

2Uij  ™*Kl  +  2  'JWij  2  *iJwi-l,J  ’ 


single  or  double  bars  connoting  intermediate  values  of  the  arguments  as 
occurring  in  Taylor's  theorem  with  remainder.  The  bracketed  part  of  the  co- 
efficient  of  w^  is  required  to  be  positive,  as  previously,  but  the 
coefficient  of  w^j  ,  unlike  before,  is  to  be  merely  non-negative.  Hence, 

0  now  is  to  be  such  that  miri  F"  -  20(1+M)  >  0  and  1  -  0(F'*+B*+C*)  >  0  , 
the  star  again  signifying  maxima  of  the  quantities  in  question. 
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Since  FN  =  1,  the  first  of  these  conditions  is 

(5)0  «  <  1/2  (1+M)  =  |  =  .166+  . 

The  second,  which  we  write  as 

(1+M)2  e2  +  (1  +  M  +  4MB)«  -  1  <  0  , 

reduces  to  (3)  for  B  =  0  and  to  the  condition 
o 

(5>i  «  <  . -  *  .156  for  B*  1/4  . 

No  stronger  condition,  as  we  shall  now  see,  arises  out  of  the  further 
argument .  Let 


(v. 


.-l,j’  WiJ»  0)  *  Nj  = 


(0,  sax  wij)  . 


From  the  above  inequality,  we  have 


<6)  “l.jtl  -  (1-'8ci3)“iJ  +  «0lJNj-f  k“iJ  - 


where  c  =  min  F"  -  2©(1+M)  =  1  -  2©(1+M).  (In  this  type  of  argument  we 
follow  Vvendenskaya.)  Since  ■©  is  supposed  to  satisfy  the  conditions  of 
the  previous  paragraph,  c  >  0  . 

Next,  like  Vvedenskaya  [6],  consider  the  quadratic  expression 


H(y)  *  H± J  (y)  =  (1  -  «C1j)y  -  |ky2  . 

H  will  be  monotonically  increasing,  i.e. ,  H'(y)  =  1  -  -  key  >  0  , 

if  y  <  (1  -  ■0Cij)Ao  .  Hence,  in  particular 
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M  H^j)  <  H(Nj)  , 


if,  for  all  i  ,  Nj  <  (1  -  ■OCij)Ac  •  In  (5.10)  we  have  proved,  however, 
that  Nj  <  sup  wiQ  <  m/h  ,  m  being  a  bound  for  v^Q  .  (In  the  present  case, 


w^j  =  W^j  as  used  on  p.  14,  and  b  =  0. )  Hence,  Nj 
condition  (7),  if  ■©  is  so  fixed  that,  for  all  i  , 
Substituting  in  this  for  c  and  also  replacing  C 


■^(1  +  M)  +  2MB  ,  we  arrive  at  the  condition 


ij 


satisfies  the  desired 
m/h  <  (l-ec^/kc  . 
by  its  upper  bound 


(l+M)(2m  -  -1(1  +  M)>G2  -  (m+2MB)«+l>  0  , 


which,  for  m  =  M  =  2  becomes 

«2  -  (2  +  4B)«  +  1  >  0 

For  B  =  0  or  B  =  l/4,  the  quadratic  expression  on  the  left  is  definite: 
the  inequality  is  satisfied  for  every  value  of  -6  and  does  not  really  con¬ 
stitute  a  restriction. 

We  are  enabled  by  (7)  to  replace  w^j  in  the  right  member  of  (6)  by 
Nj  and,  thus,  to  obtain 

’'i.j.i  s  “i  -  lkNj2  > 


and,  therefore, 


(8)  »j+1  <  Ij  -  |k»j2  . 

From  this,  it  is  easy  to  produce  an  upper  bound  for  Nj  completing  this 
alternative  proof  of  Theorem  5.2.  We  shall  show, to  be  specific,  that 
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(9) 


Nj  <  2/cjk  . 


Set  z(t)  =  2/ct  ,  z^  =  z(jk)  .  Since  dz/dt  =  -cz/2  ,  integrating 


we  obviously  have 


U+l)k 


Vl  "  Zj  =  "2 


jk 


z(t)2  dt  >  -  |kZj2  , 


or 

(10)  zj+1  >  H(Zj )  , 

where 

H(y)  5  y  -  cky2/2  . 

From  this  we  shall  prove  by  induction  that 

( 11 )  ^  Rj  »  j  =  1»2, . . .  i 

these  inequalities  are  equivalent  to  (9).  The  first  step  is  to  note  that 
z1  >  ^  ,  z^  >;  >  consequences  of  the  fact  discussed  in  a  previous 

paragraph  that  l/ck  >  m/h  >  .  Next  observe  that  H(y)  monotonically 

increases  for  y  <  l/ck  and,  furthermore,  that  z^  <  l/ck  for  j  >  2  . 

Therefore,  if  Zj  >  N^  for  some  index  j  >  2  ,  for  that  index  H(zj)>H(Nj)  . 

From  this,  (10),  and  (8),  if  for  a  particular  index  j  >  2  we  know  Zj  >  , 

we  can  conclude  that 

vi  >  u(zj)  >  »<v  >  Vi  • 

The  induction  for  (11)  is  thus  complete. 
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We  note  in  summary  that  our  first  method  bounds  4  according  to  (3) 
and  (4),  our  second  method  according  to  (3)  and  (5).  The  second  results  are 
obviously  better;  they  permit  use  of  the  values 

(12)  «  =  .166  for  B  =  0 

=  .156  for  B  =  1/4 

in  the  calculations. 

Remark:  The  left  difference  scheme  without  viscosity  of  section  2  is 
of  type  (4.4b)  with 

fij  =  Fij  *  gij  =  qij  =  sij  =  0  * 

Reasoning  like  that  above  concerning  Theorem  5.1  shows  that,  in  this  situa¬ 
tion,  M  =  m  if  «  max  F'  <  1  ,  i.e.,  if  «  <  l/(l+M)  =  l/3  .  The  further 

restriction 

(13)  «  <  l/2m  =^1/4 


is  easily  seen  to  be  required  to  bound  w^  (by  the  second  method)  according 
to  Theorem  5.2.  To  show  this,  we  start  from  (5.6)  which,  here,  leads  to  the 
inequality 


wi, j+1  -  (1  “  «Fij)wij  +  *Fj.jWi-l,j  “  2(*ijWi/ 


+  f  "  w 

+  FijVl,. 


and  thus,  in  previous  notation,  to 

.  W  2 

wi,j+l  -  wij  "  2Wij 

since  F"  =  1  and  since  the  coefficients  on  the  right  are  all  positive 
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because  of  the  first  restriction  upon  <0  .  The  previous  argument  leading  to 
the  desired  upper  bound  for  w^  will  apply  completely,  provided  we  now 
require  2m/h  <  l/k  ,  a  restriction  which  is  the  equivalent  of  (13) . 


NOL.TR  63-8 


FOOTNOTES 


—  By  generalized  solution  we  mean  essentially  a  piecewise 
continuous  solution  satisfying  the  pertinent  shock  condition 
along  discontinuities.  More  refined  concepts  of  generalized 
solution  are  to  be  found  in  [2]  and  the  accompanying  references. 

2/ 

-  For  uniqueness  theorem,  cf  [ 5l . 


3/ 

-  This  and  other  numbered  equations,  theorems,  etc.,  not  belong¬ 
ing  to  the  present  paper  are  to  be  found  in  [4]. 
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